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Abstract 

We investigate the A* (1535) —>• N'y transition in the framework of light cone QCD 
sum rules. In particular, using the most general form of the interpolating current for 
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transition form factors responsible for this channel and use them to evaluate the 
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obtained for the branching fraction is in a good consistency with the experimental 
data. 
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1 Introduction 


The investigation of the decay properties of hadrons plays essential role in understanding 
their internal structures as well as the perturbative and non-perturbative aspects of QCD. 
The theoretical study on the decay channels of the negative parity baryons, especially their 
radiative transitions, can be more useful in this regard as we have limited information 
about the decay properties of these states, experimentally. Some experimental studies on 
the photo-production and electro-production are planned to measure the electromagnetic 
form factors and multipole moments of the negative parity baryons at Jefferson Laboratory 
[1] and Mainz Microton facility [2-4], One of the main difficulty in measuring the magnetic 
moments of the excited baryons can be related to the considerably large width that they 
have. The magnetic moments of these states can be measured from the polarization observ¬ 
ables of the decay products of their excited resonances. It is hoped that the new electron 
beam facilities would allow to collect a large numbers of more precise data in the studies 
of the electro-excitations of the nucleon resonances. 

One of the main directions in obtaining essential information about the internal structure 
and natures of the negative parity baryons is to study their electromagnetic form factors and 
multipole moments both theoretically and experimentally. Such studies can also provide 
valuable information about their geometric shape. In this accordance, in the present study, 
we investigate the radiative IV*(1535) —y IVy transition in the framework of light cone 
QCD sum rules (LCSR), where IV* (1535) is the low-lying negative parity nucleon with 
J p = | . In the following, we will briefly refer to this state by N*. In particular, we 
calculate the two transition form factors ff and f p defining this channel (negative parity 
to positive parity nucleon) separating the contributions of the positive to positive and 
negative to negative parity transitions entering the physical side of correlation function. 
In the calculations, we use the most general form of the interpolating currents coupled to 
both the positive and negative parity nucleons and try to find the working region of the 
general parameter entering the general interpolating current. We also deal with real photon 
and calculate the transition form factors at zero transfered momentum squared using the 
photon distribution amplitudes (DAs). In the following, we shall refer to some studies 
on the spectroscopic and decay properties of the negative parity baryons. The mass and 
other spectroscopic properties of the negative parity nucleons have been extensively studied 
compared to their decays (for instance see [5-9] and references therein). The light-cone DAs 
of the nucleon and negative parity nucleon resonances was studied via Lattice QCD in [10]. 
The magnetic moments of the negative parity baryons have been studied via QCD sum 
rules in [11] and from effective Hamiltonian approach to QCD in [12]. The y*lV —y IV(1535) 
and 7 *N —» IV*(1520) transition form factors are studied in [13,14] in light cone QCD sum 
rules using the nucleon DAs. The electromagnetic transitions of the octet negative parity 
to octet positive parity baryons is also studied in [15] via light cone QCD sum rules. 

The layout of the paper is as follows. In next section, we derive LCSR for the transition 
electromagnetic form factors under consideration. The last section is dedicated to the 
numerical analysis of the form factors, calculation of the total decay width and branching 
fraction of the considered radiative transition. We also compare the results obtained with 
the existing experimental data. 
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2 LCSR for radiative transition of negative to posi¬ 
tive parity nucleon 

The aim of this section is to give some details of the calculations of LCSR for the transition 
form factors defining the radiative N* N 7 channel. To fulfill this aim, we consider the 
following two-point correlation function in the presence of the background photon field: 

n = i J d 4 x e ip ' x ( 0 |T (v(x) fj( 0 ) ) | 0 ) 7 , ( 1 ) 

where rj is the interpolating current coupled to both the negative and positive parity nucleon 
and T denotes the time ordering operator. According to the methodology of QCD sum 
rule, the calculation of the above mentioned correlation function is made via following two 
different ways in order to construct sum rules for the transition form factors. In the first 
way, one calculates it in terms of the hadronic degrees of freedom called as hadronic side. 
In the second way, it is calculated in terms of photon distribution amplitudes (DAs) with 
increasing twist by the help of operator product expansion (OPE) called as OPE side. These 
two sides are then matched to obtain the LCSR for the form factors under consideration. 
We apply a double Borcl transformation with respect to the momentum squared of the 
initial and final states to both sides to suppress the contributions of the higher states and 
continuum. 


2.1 Hadronic side 


The calculation of the hadronic side of the correlation function in Eq. (1) requires its 
saturation with complete sets of intermediate states having the same quantum numbers as 
the interpolating current coupling to both the negative and positive parity states. After 
performing the four-integral over x, we get 
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where • • • represents the contributions coming from the higher states and continuum. The 
matrix elements in the above equation can be parameterized in terms of residues Xn and 
A at* as well as transition form factors / 1 , / 2 , f *, /|, fj and /J at q 2 = 0 as 
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where mat and mat* are the spinors for the positive parity TV and negative parity TV* nucleons, 
respectively. The use of Eq. (3) into Eq. (2) is followed by summing over the spins of the 
TV and TV* baryons, i.e. 


22u N (.)(p,s)u N (.)(p,s) = ff + m N m . 


(4) 
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Our aim in the present study is to calculate the transition form factors f f and /J defining 
the radiative TV* —> TV 7 transition. For this purpose, we use the structures 4 4 4i i> A 4 4 
and 4 hi our calculations. The final form of the hadronic side of the correlation function is 
obtained after the application of a double Borcl transformation with respect to the initial 
and final momenta squared, viz. 

bu had = u^ AD 4 4 4+u^ AD 4 4+u^ AD 4 4+u^ AD 4+---, ( 6 ) 


where we have kept only the structures which we use to construct the LCSR for the form 
factors under consideration. The functions 11 ^ 1 ^ 3,4 are obtained as 
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where M 2 and Mf are Borcl mass parameters. As it is clear from the explicit forms of the 
functions ITf MD , they also include the magnetic moments /qv = fi + /2 and — f* + f 2 
which are calculated in [16] and [ 11 ], respectively. 


2.2 OPE Side 


The OPE side of the correlation function is calculated in terms of the QCD parameters and 
photon DAs. To this aim we need to know the explicit form of the nucleon interpolating 
current coupling to both the negative and positive parity state. It is given as 

77 = 2 e afec | (u aT Cd^j 75 u c + p (u aT C^ 5 d b ^j u c }, ( 8 ) 

where a, b and c are color indices and C is the charge conjugation operator. The parameter 
P is a general mixing parameter, which we shall find its working region in next section. 
We insert the aforementioned interpolating current into the correlation function in Eq. (1) 
and use the Wick’s theorem to contract out all the quark pairs. This leads to obtain the 
correlation function on OPE side in terms of the light quark propagators. 

The correlation function in OPE side contains three different parts: 


• perturbative part, 

• mixed part at which the photon is radiated from short distances and at least one of 
the quarks interacts with the QCD vacuum and makes a condensate, 


• non-perturbative part at which the photon is radiated at long distances. This contri¬ 
bution appears as the matrix elements ( 7 (q) | q(xi)Tq(x 2 ) \ 0 ), which are expressed in 
terms of the photon DAs having distinct twists. The T appearing in this matrix ele¬ 


ment refers to the full set of Dirac matrices, namely 
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The perturbative contributions where the photon perturbatively interacts with the 
quarks are acquired via replacing the corresponding free light quark propagator with 
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at which the free part of the light quark propagator is written as 
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where q stands for the light u or d quark. 

To obtain the non-pertnrbative contributions we replace the light quark propagator that 
interacts with the photon with 
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The remaining qnark propagators are replaced with the full quark propagators which include 
all the perturbative and non-pertnrbative parts. The full light quark propagator is given 
as [17,18] 
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where A is a scale parameter. 

As it has already mentioned, in calculations of the non-perturbative contributions, there 
appear matrix elements like ( 7 (q) | qTiq | 0). These matrix elements are expanded in terms 
of the photon DAs with increasing twists. We present the parametrization of these matrix 
elements in terms of the photon DAs together with the explicit forms of the photon wave 
functions as well as the related parameters in the Appendix. After very lengthy calculations, 
we obtain the OPE side of the correlation function in momentum space in terms of the 
photon DAs and other QCD parameters. The final form of the OPE side of the correlation 
function is obtained by applying a Borel transformation with respect to the variables p 2 
and (p + q) 2 . 

Having calculated both the hadronic and OPE sides of the correlation function in Borel 
scheme it is the time for matching the two sides in order to obtain LCSR for the transition 
form factors f 2 and /.J, defining the radiative transition under consideration. For this 
purpose, we match the coefficients of the structures /p jb /j, /p fz /q and /k from 
both sides of the correlation function. To further suppress the contributions of the higher 
states and continuum we also perform continuum subtraction. By simultaneous solving of 
four equations obtained via matching the coefficients of the above structures from both 
the hadronic and OPE sides, we get the following sum rules for the form factors under 
consideration: 
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where So is the continuum threshold and the functions nf pp (M 2 , s 0 , P, u o) are very lengthy, 
so we do not present their explicit forms here. In the above equation, m 0 and the new Borel 
parameter M 2 are obtained in terms of the Borel parameters Mf and Mf as 


M 2 

U ° ~ Mf + Mf ’ 

9 MfM 2 

M 2 = - - —— 
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Using the relation ~ we obtain Mf and Mf in terms of M 2 as 
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We also obtain uq in terms of the masses of the positive and negative parity baryons as 

m 2 N * 

uo 2 i 2 ’ 

m z N + m N , 

at which we calculate the photon wave functions and the related integrals. 

As it is clear from Eq. (13), to extract the values of the form factors ff and /.J we need 
to know the residues of both the negative and positive parity nucleons. They are calculated 
using a two-point correlation function, viz. 
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This function also is calculated via two different ways, namely via hadronic and OPE 
languages. The hadronic side of this correlation function including both the negative and 
positive parity nucleons is obtained as 
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By using the dehnitions of the residues from Eq. (3) and summing over the spins of the 
baryons N and N* and applying the Borel transformation, we obtain the final form of the 
hadronic side of the correlation function in the Borel scheme as 

BH ,had = \\ 2 N \e~ m » /M ' 2 ( $ + m N U j + |A^|e- m ^ /M ' 2 (- 0+ m N *U^), (19) 

where U is the unit matrix and M' 1 is the Borel mass parameters in the mass sum rule. 
The OPE side of the above correlation function is also calculated in deep Euclidean region 
in terms of QCD degrees of freedom and the structures and U, i.e. 

bu ,ope = nf pp (M' 2 ,4, 0 ) i> + it 2 ope (m' 2 ,4, p)u, (20) 
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where Sq is the continuum threshold in the mass sum rule. Then we match again the 
coefficients of the selected structures from both sides. By simultaneous solving of the two 
appeared equations, we get the following sum rules for the residues of N and N* nucleons 
in terms of their masses and other parameters: 


|A^| = 


e m 2 N /M' 

Tfl n + TTljy* 


|Aat*| — 


g m iv* /M' 
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Here we should also stress that the functions H'° PE (M' 2 , s' 0 , ft) and I¥° PE (M' 2 , s' 0 ,0) have 
very lengthy expressions and for this reason we do not present their explicit forms. 


3 Numerical results 

First of all we would like to obtain the behaviors of the residues of the positive and negative 
parity baryons with respect to the general parameter 0 entering the general form of the 
interpolating current discussed in the previous section. To this aim we use some input 
parameters (uu) = (d/1) = — (0.243) 3 GeV 3 , ( ss ) = 0.8 {uu) and ml = (0.8±0.2) GeV 2 [19]. 
Beside these, we need also to find the working regions of auxiliary parameters M' , Sq and 0 
entering the sum rules for the residues. The general criteria in the sum rule method is that 
the physical quantities are practically independent of the auxiliary parameters, hence, we 
look for some regions of these helping parameters such that the residues weakly depend of 
them in these regions. The continuum threshold depends on the energy of the first excited 
state with the same quantum numbers as the interpolating current for both the negative 
and positive parity states. Our calculations show that in the intervals 

1.08 GeV 2 < s ' 0 < 1.98 GeV 2 for N, 
and 

2.51 GeV 2 < s' 0 <2.73 GeV 2 for N*, (22) 

the residues weakly depend on the continuum threshold. The working region of the Borel 
mass parameter M' is obtained demanding that not only the perturbative contributions to 
the residues exceed the non-perturbative ones but also the operators with higher dimensions 
have small contributions compared to the lower dimensional operators. This leads to the 
Borel windows 

0.5 GeV 2 < M' 2 <1.5 GeV 2 for N, 
and 

2 GeV 2 < M' 2 <4 GeV 2 for N*. (23) 

It is also found that in the intervals —0.60 < cosd < —0.30 and 0.30 < cosd < 0.60 
for N as well as —0.62 < cos 9 < —0.30 for N*, at which 9 = tan~ l fl , the residues have 
relatively weak dependency on the general parameter j3. Using these working windows for 
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the auxiliary parameters as well as other inputs we find that the residue squares as functions 
of the general parameter (3 in the positive and negative parity channels can be parametrized 
as 

A 2 n (/3) = e^$i(/3), 

m AT* 

A 2 n *(/ 3) = e^ 2 (j3), (24) 

where 

$i(/3) = [2.28 - 3.20/3 + 64.10/5 2 ] x 10" 6 (GeV 6 ), 

$ 2 (j3) = [—7.38 — 4.20/3 ± 2.82/3 2 ] x 10” 5 (GeV 6 ). (25) 

Now, we perform numerical analysis of the sum rules for the transition form factors 
and use them to predict the decay width and branching ratio of the radiative N* —> IVy 
transition. To this aim, we need the photon DAs and wave functions as well as the cor¬ 
responding parameters, which are given in the Appendix. Beside these inputs and the 
residues obtained above in terms of the general parameter (3, we also need the parameters 
/ 37 = -0.00 39 GeV 2 [20] and *(1 GeV) = -4.4 GeV~ 2 [21]. 

The sum rules for the form factors /, T and /(T, contain three more auxiliary parameters: 
the Borel mass parameter M 2 , the continuum threshold So and the arbitrary parameter 
(3 entering the expression of the interpolating current. The working window for M 2 is 
obtained requiring that not only the ground state contribution exceeds the contributions 
of the higher states and continuum but the contributions of the leading twists exceed those 
of the higher twists, i.e., the series of the sum rules for the form factors converge. These 
requirements lead to the interval 2 GeV 2 < M 2 < 4 GeV 2 . We also choose the interval 
2.51 GeV 2 < So < 2.73 GeV 2 for the continuum threshold, at which the dependence of 
the form factors on this parameter are weak. Our numerical calculations also lead to the 
working regions —0.60 < cos# < —0.30 and 0.30 < cos 9 < 0.60 for the arbitrary parameter 
[3 = tan#. We depict the dependence of the form factors fl and on the Borel mass 
parameters M 2 at fixed value of so and different values of cos # in figure 1. With a quick 
glance at this figure, we see that these form factors depict weak dependence on the Borel 
parameter on its working regions. Using all input parameters and working regions of all 
auxiliary parameters as well as the masses and residues of the negative and positive parity 
baryons obtained via mass sum rules, the average values of the transition form factors fj 
and f '2 for the N* —> N 7 channel are obtained as presented in table 1. 


fl 

1.26 ±0.37 

fl 

11.18 ±3.35 


Table 1: The values of the form factors f j and /J for the N* IV 7 transition. 

At the end of this section, we would like to calculate the decay width for the radiative 
N* —> iVy transition. Considering the corresponding transition matrix elements in Eq. (3), 
we obtain the following formula for the total width of the transition under consideration: 
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Figure 1: Left: The dependence of the form factor for N* —> N 7 transition on the 
Borel mass parameter M 2 at so = 2.67 GeV 2 and different values of cos 6. Right: The 
dependence of the form factor /J for N* —> N 7 transition on the Borel mass parameter 
M 2 at so = 2.67 GeV 2 and different values of cos0. 
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Using the numerical values for the form factors and other input parameters together 
with the width of the N* state [22], we obtain the numerical values for the decay width and 
branching ratio of the radiative N* —> iVy transition as presented in table 2. In comparison, 
we also depict the existing related experimental data in this table. Looking at this table we 
see that our prediction for the branching fraction of the considered transition is in a very 
good consistency with the experimental data. 



Y(GeV) 

Br 

Present work 

(3.08 ±0.09) x 10 “ 4 

( 0.21 ±0.06) x 10- 2 

PDG[22] 

— 

(0.15 - 0.30)% 


Table 2: Decay width and branching ratio of the radiative N* —> N 7 transition. 

I 11 summary, we have calculated the transition form factors responsible for the radiative 
transition of negative to positive parity nucleon in the frame work of the LCSR using the 
most general for of the interpolating current coupling to both the negative and positive 
parity nucleons as well as the photon DAs. We found the working regions of all auxiliary 
parameters and obtained the behavior of the residues in terms of the general mixing param¬ 
eter entering the interpolating current. We used them to predict the numerical values of the 
form factors for the real photon. The values of the transition form factors are then used to 
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estimate the decay width and branching fraction of the transition under consideration. Our 
result on the branching ratio of the radiative N* —> IV 7 transition is in a good agreement 
with the existing experimental data. 
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Appendix 

The matrix elements appearing in onr calculations parametrized in terms of the photon 
DAs are written as (see [20]) : 
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(7(9)I?(®)7a*759(0)|0) = -^e q f 3l e l _ wa p£ u q a x lS due luqx ^ a {u), 
(l{q)\q(x)g s G IJ , u (vx)q(0)\0) = -ie g (qq) (£^q u - £ v q») j T>a i e l( ' a 'i +va a)i x S(a i ), 
( r r(q)\q(x)g s G tlv i^f 5 (vx)q(0)\0) = —ie q (qq) (e^q,, - £ u q^) J Va i e l Gs +va 3)^s(a i ), 
(lf(q)\q( x )g s G nu(vx) r Y a 'y5q(0)\0) = e g f 3 l q a (£^q u - £ u q^) j VoLie^ 01 ^ 01 ^* A{a.i), 
{^(q)\q(x)g s G^(vx)i^ a q(0)\0) = e q f 3 ^q a {£^q u - £ v q^) [ Pa i e i(a ’ +, ' Q!,) ' ! - T V(a ! ), 


(^(q)\q(x)a a pg s G lxl/ (vx)q(0)\0) = e g (qq) 


SX \ f 1 

£/i - q^i— [gpu - -(qpXu + q v xp) ) q D 

qx j \ qx 


£X\ ( 1 

- q,j — g av -( q a Xu + qvXa) 

qx) \ qx 


qp 


- | £u~qi 
\ 

£u - q, 


S) (' 9 “ /1 ~ qi (qaX ' 1 + %XoL 0 q0 

'fx) i 9 ^ ~ qx^^ + Qa 

£a-q a —) \g,,p - — (<2AD3 + qp x iS] g 
qx J \ qx J 

£a ~ qa qx ) ~ qx^ uXl3 + QfiXv ^) 

- (zp - (a^ - + q a xj) q» 


Vcne^+^'teTiiai) 
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Vaie^+^^iai) 


+ { £e ~ ( s ~ 


- {q V X a + q a Xv) 

qx 


% 


I 


+—(q^x„ - q v x,j)(e a qp - Epq a ) f Vaie l{a3+va3)qx %{cxi) 
qx j 

+ — (q a Zf ~ 9u*o)(£p9i, - £v%) I Da i e' < “ ,+ ”“' , ' ! “ : 7i(a i )l, 


( 27 ) 


where <^ 7 (w) is the leading twist-2 photon DAs, the photon DAs’ i/; v (u ), 0 “(m), A(aii) and 
V(cq) have twist 3; and /i 7 (m), A(u) and 71 (i = 1, 2, 3, 4) have twist 4 [20]. In the above 
relations x is the magnetic susceptibility of the light quarks. 

The measure f T>a.i entering the calculations is defined as 

/ Vai = / da q / da q / da g 5(l — a q — a q — a g ). (28) 

Jo Jo Jo 


The above mentioned photon DAs are also parametrized as [20]: 

</? 7 (m) = 6mm ^1 + 97 2 (/i)C 2 2 (u — m) j , 

■0 v (u) = 3 (3(2m- l) 2 - 1) + A (15^ - 5wf) (3 - 30(2m - l) 2 + 35(2 m - l) 4 ) , 

= (1 - (2 u - l) 2 ) (5(2 m - l) 2 - 1) ^1 + , 

A(ai) = 360a g a: 5 a: 2 ^1 + w^(7a g - 3)^ , 

V(cq) = 540w^(ag - a 5 )cK g Q; 5 Q: 2 , 

h^iu) = —10 (l + 2 k + ) C 2 (m — m), 

A (u) = 40 m 2 m 2 (3k — k + + l) 

+8(C 2 + — 3^2) [mm(2 + 13mm) 

+ 2m 3 (10 - 15m + 6m 2 ) ln(M) + 2 m 3 (10 - 15m + 6m 2 ) 1ii(m)] , 

7l (cq) 120(3^2 T C 2 )(^g CVq)c^qC^qC^gi 

T 2 {di) = 30a 2 (a 5 - aj ((m - k + ) + (Ci - Ci + )(1 - 2a 9 ) + ( 2 ( 3 -4a 3 )) , 

73(Q!i) 120(3C2 C 2 )(^<? < -^<?)5 

7 4 ( 0 ^ 1 ) = 30o: 2 (Q;q — a 9 ) ((m + ft + ) + (Ci + Ci~)(l — 2a g ) + C2(3 — 4a s )) , 

<S(c*i) — 30a g {(/i + m + )(1 — a g ) + (Ci + Ci~)(l — a g)( 1 ~~ 2a g ) 

+ C2[3(ccg - a 9 ) 2 - a s (l - a 9 )]}, 

«S(e*i) = -30a 2 {(«-/s; + )(l-a ;? ) + (Ci-Ci + )(l-a ; 7 )(l-2a 9 ) 

+ C2[3(cig - a q ) 2 - a g {l - a 9 )]}, (29) 


where, the values of the different parameters inside the wave functions are given as 972(1 GeV) 
0, uA = 3.8 ± 1.8, w* = -2.1 ± 1.0, k = 0.2, k + = 0, Ci = 0.4, C 2 = 0.3, Ci + = 0 and C 2 + = 0 
[20]. 
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